Cosmological Lower Bound on Dark Matter Masses from the Soft Gamma-ray 
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Motivated by a recent detection of 511 keV photons from the center of our Galaxy, we calculate the 
spectrum of the soft 7-ray background of the redshifted 511 keV photons from cosmological halos. 
Annihilation of dark matter particles into electron-positron pairs makes a substantial contribution 
to the 7-ray background. Mass of such dark matter particles must be < 100 MeV so that resulting 
electron-positron pairs are non-relativistic. On the other hand, we show that in order for the 
annihilation not to exceed the observed background, the dark matter mass needs to be ^ 20 MeV. 
We include the contribution from the active galactic nuclei and supernovae. The halo substructures 
may increase the lower bound to > 60 MeV. 



Recently, a narrow 511 keV 7-ray line from the cen- 
tral part of our Galaxy has been detected and mapped 
by the SPI spectrometer on the INTErnational Gamma- 
Ray Astrophysics Laboratory (INTEGRAL) sateUitef^. 
This line should be produced by annihilation of non- 
relativistic electron-positron pairs, and one of the possi- 
ble origins is the dark matter particles annihilating into 
electron-positron pairs 0, which explains the measured 
injection rate of positrons and morphology of the sig- 
nal extended over the bulge region. Therefore, 511 keV 
lines from our Galaxy as well as other galaxies may pro- 
vide a smoking-gun for the existence of light dark mat- 
ter particlesQ- Motivated by this idea, we calculate the 
spectrum of possible hard x-ray/soft 7-ray background 
(at hiy ~ 10 — 511 keV) from the dark matter annihilation 
within cosmological halos distributed over a large redshift 
range. p| and 5] have considered similar calculations at 
much higher energy (^ 1 GeV), assuming that the signal 
is coming from annihilation of neutralinos whose mass is 
typically of order 30 — 1000 GeV; thus, their calculations 
cannot be applied to the soft 7-ray background which we 
consider in this paper. Also, we improve upon the pre- 
vious calculations of this sort by explicitly including the 
free-streaming and the Jean's mass scale of dark matter 
particles, which will naturally (rather than arbitrarily) 
set the lower bound to the mass of dark matter halos. 

The soft 7-ray background has been measured by the 
A4 low energy detectors and the A4 medium energy de- 
tectors on the HEAO-1 (High Energy Astronomy Ob- 
servatories 1) satellite as well as a balloon experiment!^. 
Therefore, by comparing our calculations to observations, 
one can constrain the properties of dark matter. As the 
line is narrow, the signal must come from non-relativistic 
electron-positron pairs; thus, the mass of dark matter 
particles must satisfy 511keV < mx < 100 MeV0. We 
calculate the background intensity, I^, as 



c / dz P^,{[1 + z]!^, z) 
47r 



i7(z)(l + z)4 



(1) 



where u is an observed frequency, H{z) is the expansion 
rate at redshift z, and P^(!/, z) is the volume emissivity 



of 511 keV photons in units of energy per unit time, per 
unit frequency and per unit proper volume: 



Pi. ^ S ((1 + z)u - i^5ii) 0511/11^511 (crw)n|-. 



(2) 



Note that an electron-positron pair creates two 511 keV 
photons, and we assume that the universe is optically 
thin at hv > 10 keV. A parameter 0511 determines 
the fraction producing an electron-positron pair per one 
dark matter annihilation process. If annihilation occurs 
predominantly via positronium formation, then Q!5ii = 
1/4. While it has been shown that the positronium 
formation is the dominant process at the center of our 
Galaxy we assume that annihilation occurs directly, 
"511 = 1, for other galaxies. Here, hvsn = 511 keV 
and {(Tv) is the thermally averaged annihilation cross 
section. For flxh'^ = 0.116 j^] one finds {av) = 
[3.4,2.4,2.9] 10"26 cm3 s'^ (or [1.1, 0.80, 0.97] pb) for 
mx = [1, 10, 100] MeV, respectively (e.g. Eq. (1) in 
fioi). Here, we assume that {av) is velocity-indepedent 
(S-wave annihilation). Boehm et al.0] argue that the S- 
wave cross section overpredicts the flux from the Galactic 
center; however, we argue that it is still consistent with 
the data for mx ^ 20 MeV and p oc r~'^ * (or shallower). 
This profile explains morphology of the Galactic signal ^] , 
and one then finds that (av) ~ 0.8 pb(mx/20 MeV)^ is 
allowed at the 2-a level (see Eq. (9) in Q where the r.h.s. 



should be multiplied by 37.6/17.3 ~ 2.2 for p (x 



As we show later, constraints from 7-ray background fa- 
vor this parameter range, and our results do not violate 
the Galactic constraint. 

Since the number density of dark matter particles, nx, 
is usually unknown, we use the mass density, px = 
nx Tnx , instead of nx ■ After multiplying hy = 
J^5ii/(1 + z), one obtains 
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where 1 + z — v^u/v, {Px) ^ is Pjc averaged over 
proper volume at z, and Cx{z) = {p^x) ^ I {px)\ is the 
dark matter clumping factor. (We have used {px) ^ — 
10.54 Ox/i^(l + keV cm"^.) It is worth noting that 
since {av) oc (fijc /i^)"-'^ j^, the intensity scales approxi- 
mately as vly oc \Cx{z)/m\'\ {Ux /^ll^^ /i(l + z)^/^ for 
z > 1; thus, the larger the mass per dark matter particle 
is, the smaller the predicted 7-ray background becomes 
(oc m'^), and cosmology and redshift (or v) dependence 
of the signal is almost entirely determined by Cx{z), i.e., 
how dumpiness of the dark matter halos evolves with 
redshift. Throughout this paper, we adopt the best-fit 
power-law ACDM model from the first year Wilkinson 
Microwave Anisotropy Probe {WMAP) satellite ,9]. 

Eq. allows one to calculate the intensity directly 
from A^-body simulations which, in principle, give the 
full evolution of Cx{z)', however, such simulations are 
computationally challenging as they must resolve indi- 
vidual halo profiles at small distances (e.g., ^ kpc) in 
a large simulation box (e.g., 1 Gpc), which requires a 
large number of particles. Instead, we take an alternative 
approach. The clumping factor is essentially determined 
by the clumping of individual halos and the volume occu- 
pied by halos (or the number of halos). The former can 
be obtained by simulating individual halos using high- 
resolution simulations with smaller box sizes, while the 
latter can be obtained by simulating many halos using 
low-resolution simulations with larger box sizes. There- 
fore, by combining these simulations, one can calculate 
the clumping factor. This somewhat empirical approach, 
called a halo approach (see 12J for a review), turns out 
to be very powerful in making accurate predictions for 
clustering of dark matter particles in a highly non-linear 
regime. By using this approach, we calculate Cx{z) as 



Cx{z) = 



{l + zf 
{Px)l 



dM 



(4) 

where px{M,r) is a dark matter halo profile, and 
dn{M,z)/dM is the comoving number density of ha- 
los in the mass range of M ~ M + dM at z. For 
dn/dM we adopt an empirical fit to A^-body simula- 
tions obtained by with a — 0.75, p = 0.3, and 
A = 0.322 (see [l^ for the detail). Why does the integral 
in Eq. (0} have a minimum mass? When the length scale 
is smaller than the free-streaming scale and/or Jean's 
scale of the dark matter particles, density fluctuations 
are suppressed and halos do not form. Although the 
actual effect is not an abrupt cut-off and will cause a 
gradual suppression of halos toward low mass scales, we 
approximate it as the low mass cut-off in the integral; 
thus, no halo would form below Mmin = max{Mp, Mj), 
where Mp is the free-streaming mass and Mj is the 
dark matter Jean's mass. For Mp, we adopt the expres- 
sion obtained by Q: Mp « {'iTr/3)'p:;^{\p/2)^ , where 



\p = 2.72(f7x/i^)^/^(mx/lkeV)-4/3 Mpc, and p^;: is the 
present-day mean matter density. With this definition 
one obtains Q 



Mp ^ 8.3 • 10"^ Mq 



0.13 / V 0.11 / VlMeV 



The dark matter Jean's mass is given by [lE 



(5) 



3/2 

Mj = 38.79 Mq{^) {^) [ ^ ) , (6) 



mx /Xi?\3/2 / I + z 
Mev) vT2/ V 3069 



for 1 + z > 1 -|- Zeq ~ 3069, which is a valid expres- 
sion for cold thermal relics which decouple from ra- 
diation at xp — mx/Tp > 3 {xp « [11 — 15] for 
mx 100] MeV; e.g. see [l3). 

Eq. (0J can be cast into the product of the collapse 
fraction and the mean "halo clumping factor" : 



Cx(z)- A.Feon(z). [C^-i°] 



(7) 



where A is the mean halo overdensity of a halo 
in units of the cosmic mean matter density, the 
collapse fraction Fcon{z) = / dAIj^M/'p^ is the 
mass fraction collapsed into cosmological halos at 
z, [A] = J dM-^MA/ f dM-^M, and C^^'° = 

/ d'^r ^^^^ — ^ / / d^r is the "halo clumping factor" 

defined in terms of the halo mean density {px)haio ^ 
J d^rpx/ J d^r ^ A ■ {px)^- As illustrated in Fig. |(2J), 
the early time Cxiz) is mainly determined by Fcoii{z), 
while the late time Cx{z) by [C^^'°]. As it can be seen 
from Eq. Cx{z) is very sensitive to the density pro- 
file, whose properties are not fully understood (or ob- 
served) yet. In order to quantify uncertainties associated 
with the density profile, therefore, we adopt two different 
models for density profiles. 

Case A: The Navarro- Frenk- White (NEW) Profile. 
The NEW profile is an empirical fit to radial profiles 
of dark matter halos in A^-body simulations ,16i] . This 
profile has a central cusp, px oc r~^, and is therefore ex- 
pected to produce large annihilation signals. The NFW 
profile is given by px{t) = Ps(r/rs)~^(l -I- r/rs)"^^^. 
The mass M enclosed within the virial radius, rvir, is 
M — 4:TTpsr^ [ln(l + c) — c/(l + c)], and the scale radius, 
Ts, is rs = J'vir/c, where c is called the concentration 
parameter. With these definitions, one obtains 



Chalo 
X 



C3(l_l/(l + c)3) 

9(ln(l + c)-c/(l + c))2 



(8) 



For instance, C^^i°(c = 3) = 7.3, C^^^°(c = 10) = 50, 
C^^°{c = 20) = 203. For A, we use an approximate form 
A = (187r2 + 82x - 39x^)/n{z), where n{z) is the ratio 
of mean matter density to the critical density at z, and 
X — 0(z) — 1. This expression is valid for a flat ACDM 
universe [l^lisl ]. The concentration of dark matter halos 
found in A^-body simulations has a log-normal distribu- 
tion with a median value of |l8j ] 



c(M, z) 



l + Zc 

' l + z 



(9) 
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FIG. 1: Predicted spectrum of the soft 7-ray background 
from cosmological halos with the NFW profile. From the 
top to the bottom, each set of curves represent mx = 
1, 10, 100 MeV. The solid lines show the canonical concen- 
tration model [Eq. @], while the dotted, short-dashed, long- 
dashed and dot-dashed curves show the models with upper 
cut-offs of Cmax =50, 20, 10 and 3, respectively. Signals with- 
out upper cut-offs are well fit by a power-law for v > 40keV, 
as indicated by the uppermost line. 



where the collapse redshift, Zc, is implicitly given by a re- 
lation Mf,{zc) = 10^^ M. {M^,{z) is the non-linear mass 
at z.) As lower mass objects collapse at higher Zc, the 
concentration decreases as M. (The lower mass objects 
have steeper profiles.) Following (l^, we take into ac- 
count a log-normal distribution of c with the dispersion 
of a{\n{c)) ~ 0.2. While we use Eq. for all the range 
of M and z in our analysis, we should keep in mind that 
this fitting formula is valid only for a limited range of M 
and z covered by iV-body simulations. 

Fig. n shows the predicted soft 7-ray background for 
the NFW profile. The predicted signal, I'l^, is roughly 
proportional to i^^-^, which simply reflects the fact that 
the universe becomes more clumpy at lower z as C{z) oc 
(1 -I- z)-^-^ (recafl that hi/ = 511 keV/(l -|- z)), as seen 
in Fig. 12 One should, however, keep in mind that we 
have simply extrapolated the fitting formula for c(M, z) 
[Eq. ((5J)] to the regime where simulations are no longer 
valid. For the mx — 1 MeV case, for example, the 
free-streaming mass gives the maximum concentration 
of c ~ 70 at z = 0. Since we don't fully understand 
halo profiles at such low mass or high concentration, we 
also apply an arbitrary, hypothetical upper limit to c 
and investigate sensitivity of our results to the change in 
concentration parameters. This toy model will substan- 
tially reduce the contribution from small mass halos with 
c{M,z) > Cmax- As c{M , z) OC (1 + z)^^ , the signal at 
lower z (higher hv) is suppressed. As a result, i^/i, be- 
comes almost flat for v > v^:, where corresponds to a 
redshift, z*, after which the clumping factor stops evolv- 
ing fast and evolves only slowly as (1 + z)~^/^ or even 
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FIG. 2: Upper panel: The clumping factor Cx (z) for different 
mx's. Solid lines, from the bottom to the top, represent 
the clumping factor for mx ~ [1, 10, lOOjMeV respectively. 
Dotted lines represent the clumping factor for Cmax ~ 3, also 
with the same meaning from the bottom to the top. Two 
power laws reflecting the asymptotic behavior of Cx{z) are 
also plotted. The trend for different mass curves indicates the 
fact that the minimum halo mass is smaller for larger mx, 
thus obtaining higher clumping. Lower panel: The mean halo 
clumping factor [Cx'''°(2)l and the collapse fraction -Fcoii(z) 
for mx = 1 MeV. The solid line is for no upper limit on c, 
while the dotted line is for Cmax = 3. It clearly shows that 
the low redshift behavior is determined mainly by [Cx^'°(z;)], 
while the high redshift behavior by _Fcoii(z). 



slower (Fig. ^ ^ . Sensitivity of the predicted spectrum 
to the concentration parameter model is demonstrated 
more in Fig. |51 For the canonical concentration parame- 
ter model given by Eq without any upper limit, the 
smallest mass halos always dominate; on the other hand, 
once the upper limit on c is imposed, the largest con- 
tribution comes from M ~ M*(z), effectively removing 
contribution from the lower mass halos as the lower mass 
halos hit Cmax earlier than the higher mass halos. 

Case B: The Truncated Isothermal Sphere (TIS) Pro- 
file. The TIS 20] is an analytical model of the post 
collapse equilibrium structure of halos resulting from the 
collapse of a top-hat density perturbation. The TIS is 
the minimum energy solution to the Lane-Emdcn equa- 
tion, matching energy of the TIS to the energy of the 
top-hat perturbation. The TIS density profile has a soft 
core and thus provides a much smaller annihilation sig- 
nal than the NFW profile. Note that the TIS fits the 
observed rotation curves of dwarf spheroidals and low 
surface brightness galaxies fairly well, as opposed to the 
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FIG. 3: Fractional contribution from different NFW halo 
masses at different redshifts. The left panel shows the canon- 
ical concentration model [Eq. Q] with no upper limit, while 
the right panel shows the model with Cmax = 3. In both pan- 
els, the solid, dotted and dashed lines correspond to z —0, 4 
and 50 respectively. 



NFW. The TIS density profile is given by |2C 



21.38 



19.81 



9.082 _j_ (r/ro) 14.622 + (r/ro) 



for r < rt- Here, Pm{z) is the mean matter density at z, 
Tt is the truncation radius, 



rt 



M 



0.187 



l + z\ W^Mr. 



(a^/l2)-l/3 Mpc, (11) 



and ro is the core radius given by tq = r(/29.4. We find 
that the predicted signal for the TIS model is very similar 
to that for the NFW with Cmax = 3, which is consistent 
with the TIS profile being much flatter (i.e., much less 
clumpy) than the NFW with the canonical concentration. 
This can be more easily explained by Eq. [7| At z ^ 1 , 
A « 178 and Cx{z) « 178 x 7.3 x Fcon for the NFW with 
c^^^ = 3, while the TIS has A w 130, C^'^'" « 11.5, and 
so Cx{z) « 130 X 11.5 X Fcoii- 

AGNs and Supemovae. There are other sources for 
the soft 7-ray background. The Active Galactic Nuclei 
(AGNs) are probably the most dominant contribution 
up to 100 keV, and it may continue to dominate up to 
~ 500 keV|2ll l2a|. depending on a cut-off energy scale 
of the AGN spectrum, Ecut, which is a free parameter. 
The other candidate source is the Type la Supernovae 
(SNIa). 'E] have calculated a spectrum of the soft 7-ray 
background from SNIa, finding that the SNIa contribu- 
tion can account for the observed intensity. Our recent 
calculations, however, show that they overestimated the 
signal by a factor of ~ 10 (they assumed ~ 10 times larger 
supernova rate than observed). The SNIa contribution 
is thus negligible at the energy scale of our interest 0. 
Figure ^ compares the dark matter annihilation, the 
Compton-thin AGNs with Ecnt = 500 keVjQ, and the 
supernovae contribution |2^, as well as the observed soft 
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FIG. 4: Comparison of the dark matter annihilation (dashed 
lines for the NFW and dotted lines for the TIS), the 
Compton-thin AGNs 22] (solid), and the type la supernovae 
contribution j2^j (dot-dashed), as well as the observed soft 7- 
ray background|g (points with errors). The predicted total 
signals are plotted as the thick dashed and dotted lines, which 
set upper limits on the dark matter mass from the measured 
points. 



7-ray background foj. We find that an acceptable range 
for the dark matter mass is 20 ^ mx ^ 100 MeV for the 
NFW profile. (Note that the upper limit comes from re- 
quiring that electron-positron pairs be non-relativistic.) 
The uncertainty includes statistical uncertainty as well 
as potential systematic uncertainties associated with in- 
strumental calibrations. It has been shown that the back- 
ground signals measured by various experiments show a 
large scatter (e.g., "25]), and it may be possible that the 
measured background is uncertain up to 30% 22] . 

One may relax this constraint in various ways. If 
the AGN spectrum has a smaller energy cut-off, i^cut ~ 
100 keV, the AGNs do not contribute to the 7-ray back- 
ground for hi/ > 100 keV, and mx ^ 12 MeV becomes 
acceptable (Fig. \^ . Bigger uncertainties come from the 
halo model. If we adopt the TIS or the NFW with 
Cmax — 3, we obtain mx ^ 5 and 1.7 MeV for with and 
without the AGNs contribution, respectively (see Fig. 01 
and[SJ). Such a low concentration parameter or filtering of 
low mass halos may be possible if dark matter has some 
finite self-interacting cross section[2^l2^. In addition to 
the average shape of px , the clumpy substructure within 
halos would contribute to the signal sig nificantlyJlll^. 
While the abundance and properties of substructure are 
uncertain, the analysis in P| shows that the substruc- 
tures can enhance the clumping factor by more than a 
factor of 10. As the intensity is proportional to m~^^, 
the contribution from the substructures can increase the 
lower limit on mx by more than a factor of 3, giving 
mx > 60 MeV for the NFW profile. 

In summary, we have calculated the soft 7-ray back- 
ground of the redshifted 511 keV photons from the dark 
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which is consistent with that from the Galactic center Q 
for p oc r'~°'* or shallower profile). Recently, Bea- 
com, Bell and Bertonel^^ have shown that the internal 
bremsstrahlung of electrons and positrons contributes to 
the 7-ray background at > 1 MeV. Using COMPTEL 
and EGRET data, they have obtained mx < 20 MeV. 
Putting all the constraints together, mx ^ 20 MeV may 
be favored. However, in order to have a more robust 
constraint, further understanding of cosmological non- 
linear structures is required. Obtaining a more reliable 
model for AGNs is equally important. Uncertainty in the 
cut-off energy, Ecut, is the dominant uncertainty in the 
predicted soft 7-ray background signal at ^ 100 keV. 
If the cut-off energy is smaller than ~ 100 keV, the soft 
7-ray background signal may serve as evidence for the ex- 
istence of dark matter particles with mx ~ C(10) MeV. 



FIG. 5: Dark matter annihilation only. Smaller dark matter 
masses are allowed, (see also Fig.^]) 



matter particles annihilating into electron-positron pairs 
in cosmological halos. Our fiducial model based on the 
current iV-body simulations, compared to the observed 
7-ray spectrum, places limits on the dark matter mass, 
20 ^ mx ^ 100 MeV. (If the cross section is allowed 
to be free, we find {{av) /ph){MeV /mx)"^ < 2.1 x lO'^, 
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